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Abstract 


This work describes the fundamental mechanism arising during the bursting of a bubble at 
an air-liquid interface. A single bubble is followed from an arbitrary depth in the liquid, 
up to the creation and motion of the film and jet drops. Several phenomena are involved and 
their relative order of magnitude is compared in order to point out the dimensionless para- 
meters which govern each step of the motion. Furthermore, this study is completed by high 
speed cinematography. The characteristic bubble radius which separates the creation of jet 
drops from cap bursting without jet drops is shown to be a = (c/pg)*/2, The corresponding 
numerical value for water is 3 mm and agrees with experimental observations. 

Introduction : physical mechanisms description 

The burst of gas bubbles on a liquid interface leads to a mass transfer in stratified 
two-phase flows which is of great importance in many fields such as ocean-atmosphere ex- 
changes, aerosol generation, boilers, dega 2 eif ication processes, etc... The physical me- 
caniims of such a phenomenon are still not wellknown. The burst occurs during a few tens of 
microseconds and involves many physico-chemical effects. The film cap draining is due to 
gravity and surface tension and its thickness (in the range of 0.5 to 10 -: micrometers) may 
be small enough to produce non negligeable intermolecular forces which depend on the impuri- 
ties contained in the liquid phase. The film cap breakage generates an aerosol which is 
usually split into two distinct families which are the "film drops" and the "jet drops". 

. The film drops are created by the tearing of the interfacial cap which propagates frcm an ini- 
initial hole at a very high speed (8 m/s) . Due to the tangential motion of the liquid film, 
some drops are ejected horizontally and the others are blown vertically by the air which was 
initially inside the bubble. Another possible mecanism is the sudden shattering of the whole 
cap which bursts into droplets in any direction . Film drops are very small (1-20 pm) 
and their number increases with the initial bubble diameter 0b (about lOO when 0b = 2 mm, 
and lOOO when 0 = 6 nun) . They are ejected up to a height in the order of 10 to 20 mm above 
the free surface level. 

. The jet drops are produced by the very fast motion of the bubble bottom interface which 
is in a state of non-equilibrium when the film cap breaks. Instantaneously, gravity and sur- 
face tension act to create a strong interface deformation and give rise to the formation of 
a jet in the upward direction. Because of instabilities, the liquid jet is fragmented into 
several drops (from 1 to 5 depending on the bubble diameter 0b, in the range 6 mm to 0.1mm). 
For values of 0^ less than 2 mm, the ejection height of jet drops may reach 100 0b above the 
free surface level. 

These two families lead to a transfer in the gas of two different kinds of chemical sub- 
stances, namely the surfactants spread on the free surface (film drops) and the particules 
j.n suspension or dissolved inside the liquid (jet drops). 

The aim of this paper is to consider a dimensional analysis of aerosol production in or- 
der to reach a better understanding of the phenomena involved during the bursting of the 
bubble. The notion of dimensionless parameters, characteristic times and predominant mecha- 
nisms may be used to define simplified experimental studies and to search for approximate 
analytical solutions. In what follows, several effects are analyzed : 

- underwater bubble rising - film drainage - cap tear propagation 

- film drops ejection - free surface motion after film cap rupture. 

Underwater bubble rising 
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The dimensionless form of equations governing the relative liquid flow past the bubble 
is written with the characteristic bubble size a, the pressure variation pag (involving 
the liquid density p and the acceleration g gravity) and an unknown characteristic time t 
which leads to the reference velocity a/t. The bubble motion is obtained from the fundamen- 
tal dynamic law : 
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The bubble shape is deduced from the normal stress equilibrium 
P a 
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n is the unit normal vector, p., C, v and o respectively stand for the internal bubble pres- 
sure, the mean curvature, the Aynamic viscosity and the surface tension coefficient. 

The local pressure is defined by the relation 

P a 

P = + z^ft) + P, 

* pga B ^loc. 

in which, for the most part, the constant p /pga is much greater than unity. 

cl 

It is easily seen, from equation (2) . that high Reynolds number flows correspond to the 
characteristic time evolution ig= ( a/g) */ 2 , and situations dominated by viscous effects 
are scaled by the reference time t v = v/ag. The relative Reynolds number : 


= a g/v = (a/aj) 


(4) 


leads to a characteristic length a 3 = (v 2 /g) 1//2 [^50 ftm for water). In equation (3), ano- 
ther length ai = 2(o/pg) 1 / 2 (^6 mm for water) is defined from the dimension) ess number : 


l = 4o/pga 2 = (a 2 /a) 2 = <T g /t s ) 2 (5) 

The new characteristic time t_ = pa 2 /4o concerns the interface vibrations due to the surface 
tension. Depending on the bubble diameter a, the following situations occur. 

. a >> a 2 , R r >> 1, motions due to surface tension are frozen, but the bubble can oscil- 
late with instationnary terms of equation (3) (bubble wake for instance) on the time scale 

T g‘ 

. a! << a << a 2 , R r 1 but *"he non stationnary terms in equation (3) are negligeable 
and the bubble keeps its spherical shape. 


. a << aj, Rr << 1 and the comparison between the last two terms of equation (3) leads to 
a new length a 3 = pv 2 /2o ('vS.lO - ' u m for water). Bubble diameters less than a 3 will not be 
considered in this paper. Therefor surface tension is the dominant phenomenon and deforma- 
tions cannot occur. We will write * * a 3 /a * pv 2 /2«a (always less than unity for water) . 


The characteristic velocities of the rising motion are 
(stokes drag) when a < a*. 
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Film draining 

There are two fundamental mecanisms which cause the film draining, when the bubble ar- 
rives on the 'ree surface. First, the gravity which creates the bubble motion and inertia 
terms (when a > a 3 ) increases the local pressure p a on the free surface. On the relative 
stagnation point, this increase is of the order of pga. Secondly, the free surface curva- 
ture leads to a pressure increment in the order of 4o/a. Comparison between both phenomena 
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is given by the value of I » 4o/pa 2 introduced in the last paragraph. Consequently, for 
bubble diameters larger than a 2 , the effect of gravity is predominant, and for aj << a << a t 
the surface tension plays the most important role. 

A very crude analysis for z << l shows that (from the continuity equations) the liquid 

pushed up by the rising of the bubble (order £ a 2 /ag) must be evacuated through the lateral surfa- 
ce of a cylinder (it a h /ag) . This balance gives an order of magnitude of the height below 
the free surface which must be reached by the bubble top before undergoing a non negligeable 
free surface deformation. When a^ << a << , the characteristic bubble rising velocity /ag 

has not the same order of magnitude as the characterirMc velocity due to the surface ten- 
sion pressure drop : (4a/pa)*'* and the same crude analysis leads to : 

j a 2 /ag = 2 i a h (o/pa ) 1/2 and h * — ^ 7 ? << a 
For very small bubble '^3 << a << a^), the same trivial argument gives : 



The characteristic 1 aining time is easily deduced from the above considerations, and 
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from the volume to be drained (■v ~g- when 1 << 1 ; - - when a^<< a << a 2 ; 


and 
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. As an example in each case, one obtains respectively for a = 1 cm, 1 mm and 


10 um, the following values : = 0.07 s, 4 ms and 0.5 ms, in water. The end of drainage 

in the first two cases is governed by viscous effects when the Reynolds number based on the 

1/2 1 /2 

the film thickness e is of order unity i.e e * v/(ag) when z << 1 and e v(pa/o) 
when z >> 1 . These limiting thicknesses in water are respectively 3 pm and 0.3 pm when 
a = ’ cm and a = 1 mm. At this stage, it is necessary to be very careful about the boundary 
conditions to be written on the film boundaries. For very dirty liquids, Haberman and 
Morton have shown that a vanishing velocity condition must be considered. In any case, one 
must take into account the influence of Van der Waals intermolccular forces which are, for a 
clean liquid, given under the form of a non isotropic tension in the normal direction to the 
film. Its value for a thickness e is (Sheludko) : 

x = — ‘yP-- ‘ (MKSA units! 


this influence is non-negligible when the pressure gradient is of the same order : 

4 10*" 1 5 10 "^ 

a > a 2 ( F < * 1 ), eg ^ — : — 2 and e ^ — : — \TT = 0 . 5 um 

e (rg) ' 

a * a- ( i •>> 1) , ~ - 3 4 - 1P .-1 and e ^ 10" ? a 1/2 M.K.S.A. (0.3 pm when a - 1 mm) . 

‘ a^ r l/ * e 4 

These orders of magnitude for the film thickness is in agreement with the experimental eva- 
luations. 
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Cap tear propagation 



Assuming that a hole is created instantaneously, the driving force is due to surface tension 
and leads to a characteristic velocity : v f ■ (a/pe) 1 / 2 . If the characteristic thickness 
e 1 pm is admitted, it follows that in water (o ^ 0.0'<) the tear propagates with a veloci- 
ty of 8 m/s which is exactly the value observed experimentaly with a high speed camera. This 
evaluation ii valid, only when the liquid flow inside the film is dominated by the inertia 
terms. This happens for a local Reynolds number greater than unity : 
v,e 2 

R, = — >> 1 -> e >> — ( = 0.01 pm) 

IV 0 


So, it is obvious that the experimental study of cap tear breaking may be performed without 
any precaution because any other phenuiienon is frozen with respect to cap tear propagation. This 
has been done in the movie presented at the end of this paper. 

This experiment shows the cap tear of a large bubble (0^ = 4 cm) at rest on the free sur- 
face of ordinary drinkable water. The hole is created on the circle which delimits the bub- 
ble cap on the free surface. It is due to a much smaller bubble sticking to the larger one 
and which causes a point of weakness. Soma papers (Vrij , Ivanov ,...) explain the film 
breaking by the non stable modes of capillarity waves, but this is probably true for bi- 
distilled water only. 


Film drops ejection 

The very complex nature of film cap fragmentation is not investigated in this paper, but 
we will consider the behaviour of a film droplet assumed to have a size in the order of 1 pm 


The blowing-out is due to the fact that the pressure level Inside the bubble before brea- 
king is larger than on the outside. The sudden pressure drop is of order o L ag when I << 1 
and o/ a when l >■> 1. For high Reynolds number gas flows, this produces velocities of order 
(PLag/.^) 1 ' 2 or (a/a ) l ' 2 . For each case, the gas Reynolds number is : 


(f7) 


1/2 


a 3/2 

>> 1 if a » 5 pm 

V G 


>> 1 if a >>0.1 A 


(air an^ water) 


(air and water) 


2 

The characteristic time of a particule transport due to the viscous drag is t - r /v G f(. 
where r denotes the particule radius. For a micronic droplet this time is 1 mill - second , as 
a result, film drops follow instantaneously the air flow produced by the breaking. For bubble 
diameters : 1 cm, 1 mm, 10 pm, the initial velocity of droplets are respectively : 10 m/s, 

10 m/s and 100 m/s. 


Experimentally, it is know that a 10 pm bubble does not produce any film drops, but one 
can see that for most bubbles, there is a strong interaction between the cap tear propaga- 
tion (8 m/s) and the air ejection (10 m/s) . If the air flow is roughly supposed to be simi- 
lar to a source flow, the velocity V r at a distance r is given by V r * Vi(a/2r) . A numeri- 
cal value of V- with a = 1 mm and r * 1 cm is 2.5 cm/s (less than 1 km/h). This argument 
proves that film drops are immediatly convected by the wind in ocean-atmosphere exchanges. 

Free surface motion afeer film cap rupture 

It has been shown that film tear propagation occurs instantaneously with resr p to the 
bubble rising velocity. This result confirmed by Mac Intyre's visualisations a by the pho- 
tographs given in the present paper. As a result, the curvature of the bubble is greater at 
the bottom than the sides. This causes stronger surface tension forces on the deepest part 
of the bubble interface and produces the jet when the sur f ace tension phenomenon is not 
negligible with respect to the gravity effect. 
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Sb 


The dimensionless equations governing the free surface motion are written assuming a 
potential flow (v = 70) for the liquid phase. It will be proved ( later on, that inertia 
dominates the viscous terms. In a cylindrical coordinates system, the bubble diameter a is 
chosen as the reference length in the radial and vertical directions. The characteristic 
time T needed for the free surfa-e z = i(r,t) to reach an altitude of order a is unknown at 
this stage. 
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The reference potential and the time T are obtained from the free surface boudary condi- 
tions. Equation (6) leads in any case to the relation = a 2 /T. From equation (5), the cha- 
racteristic time T depends on the comparison between = pga 2 /4 (gravity) and the mean 

dimentionless curvature C which is of order unity, z -•> 1 : surface tension is the dominant 
phenomenon and T = ,ai 0 /'' = pa 2 /oT. The characteristic time and velocity are respectively s 

T * (pa 3 /o) 1/2 and V = (o/pa) 1/2 

i << 1 : the free surface motion is governed by the gravity effect. One obtains respectively 
T * (a/g) 1/2 and V = (ag) 1/2 


These results are valid only for high Reynolds number flows. When r >> 1 (a << a 2 , 0.6 cm 
for water) R * aV/ . = (a o/ .v)l/ 2 which is larger than unity for bubble diameters larger 
than a = pv 2 /c (0.01 m for water). When i: << 1, R » (a-^g/u 2 ) i/ 2 and the bubble diameter 
must be larger than aj ■* (u 2 /g) 1 / 2 (50 yir. for water). This condition is always fulfilled 
because a > > 0.6 cm. 




The characteristic bubble diameter a 2 ■ (4 'pg) 1 ^ 2 splits the possible free surface beha- 
viours into two classes : 
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a > a 2 : the film cap rupture leads only to free oscillations of the free surface without 
jet, but many film droplets are blown into the atmosphere. 

a < &2 '• the predominance of the surface tension effect leads to a strong driving force 
at the bottom of the bubble which produces a jet. Surface tension forces still act on the 
jet surface and cause pressure gradients along the jet axis, oriented from the smallest sec- 
tions to the largest ones. The drairing of the lowest sections due to this mechanism leads to 
the jet breaking into one or :everal drops. 



The dimensionless equation (5) shows that the gravity effect becomes important for a di- 
mensionless z of order l . So the jet top can reach an altitude of order ar. The experimental 
evaluation : lOO a, given in the introduction corresponds to a ■* 0.6 mm and an altitude of 
60 cm. Dissipations due to viscous effects are given by the Reynolds number R ■ (a a/ P v 2 )l/2 
^ 8.10 for a 0.6 mm bubble diameter. In that case, the boundary layer on the free surface, 
in which viscous draining occurs has a thickness in the order of a hundredth of a bubble 
diameter. 


Flow visualisation 


The photographs given in this section have been obtained with a high speed camera HYCAM 
(7500| images per second in the present case). The first sequence represents a cap tear which 
propagates on a 4 cm bubble diameter, from a smsll bubble K 1 mm) which was initially atta- 
ched to the large one on the free surface. One can observe that some droplets tangential- 
ly ejected by the tear motion have a much larger size than 1 micrometer. 

The next two sequences show jet drops due to a 1 mm bubble diameter. In ’JET N°l", the 
hole created instantaneously by the breaking of the fiTm is shown. The jet thickness is of 
order 0. 1 mm corresponding to velocities of order 1 m/s. The half-sphere in the up- 

per part of each view is a liquid drop sticking on the end of a 1 mm cylindrical wire placed 
there to give an indication of the length scale. The "JET N # 2 M gives more details on the jet 
breaking mechanism. One can observe that che top drop has no motion during its growth, but 
its upward motion continues after its creation. 

Concl idlng remark s 

This paper is a first step in the study of aerosols created by bursting bubbles on an air- 
liquid interface. Jet drops have been studied and visualized in the past ; however a theore- 
tical model for the prediction of tha size and number of the drops would be useful for a bet- 
ter understanding of the jet breaking mechanisms. 

Neverthless, the lack of knowledge is more important for film droplets and a lot of expe- 
rimem.il work will be necessary in order to visualize the film breaking and to predict the 
full size distribution of aerosols in a real industrial statified two-phase flow figuration. 


This study has been supported by a CNRS* contract. 


* Centre National de la Recherche Scientifique (France) . 
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